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We introduce a three-dimensional Plebanski action for the gauge group SO (4). In this 
model, the B field satisfies quadratic simplicity constraints similar to that of the four- 
dimensional Plebanski theory, but with the difference that the B field is now a one-form. 
We exhibit a natural notion of "simple one-form" , and identify a gravitational sector, a 
topological sector and a degenerate sector in the space of solutions to the simplicity con- 
straints. Classically, in the gravitational sector, the action is shown to be equivalent to that 
of three-dimensional first order Riemannian gravity. This enables us to perform the com- 
plete spin foam quantization of the theory once the simplicity constraints are solved at the 
classical level, and to compare this result with the various models that have been proposed 
for the implementation of the constraints after quantization. In particular, we impose the 
simplicity constraints following the prescriptions of the so-called BC and EPRL models. We 
observe that the BC prescription cannot lead to the proper vertex amplitude. The EPRL 
prescription allows to recover the expected result when, in this three-dimensional model, it 
is supplemented with additional secondary second class constraints. 



Introduction 

Spin foam models [IHl] constitute an exciting proposal for the definition of a background inde- 
pendent and non-perturbative quantization of general relativity. They were introduced originally 
[3] as a way to implement the dynamics of loop quantum gravity [HI E] , and can be thought of as 
representing a sum over histories of the gravitational field [8J. 

In three spacetime dimensions, pure gravity being a topological theory (in the sense that it has no 
local degrees of freedom), spin foam models can be obtained as an exact simplicial path integral for 
the first-order Palatini action, where the gauge group is taken to be SU(2) in Riemannian signature 
(the whole symmetry group is larger than SU(2) and depends on the sign of the cosmological 
constant A). This corresponds to the celebrated Ponzano-Regge |9j and Turaev-Viro [lO] models 
when A = and A > respectively. Their vertex amplitudes are simply given by SU(2) and 
Uq(su(2)) (for q a root of unity) 6j coefficients. It has been shown [11] that the Ponzano-Regge 
amplitudes can be obtained from a canonical quantization of 2 -|- 1 gravity in the spirit of loop 
quantum gravity, in which the kinematical states are given by spin networks, and the dynamics 
leads to a spin foam evolution. The Turaev-Viro amplitudes can also be obtained as the scalar 
product between physical states which are in turn given in terms of unitary representations of some 
quantum groups closely related to Uq(su(2)) [121 [13]. Therefore, in three dimensions, there is a 
clear understanding of the relationship between the so-called canonical and covariant approaches 
to quantum gravity (at least when the signature is Riemannian). 

In four spacetime dimensions, the construction of such a correspondence between the canonical 
and spin foam quantizations is still an open problem |14H16j . Since it is technically too involved 
to follow the construction of the three-dimensional theory and perform the straight simplicial path 
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integral quantization of the four-dimensional Palatini action (see [T7] for an attempt to follow this 
direction), spin foam models have to be derived using an alternative strategy. More precisely, the 
approach to four-dimensional spin foam models is based on the fact that gravity can be formulated 
as a constrained topological field theory, by virtue of the so-called Plebanski action [TH] . This latter 
is defined as the sum of a well-known topological BF action |19j , plus a set of constraints on the B 
field ensuring on-shell the condition that it comes from the exterior product of two tetrad one-form 
fields^. With this solution for the B field, the Plebanski action reduces at the classical level to 
the usual Hilbert-Palatini action of first order general relativity. To derive a spin foam model, the 
following general strategy is adopted. First, the topological BF theory is discretized exactly on a 
two-complex, and the basic variables, i.e. the B field and holonomies of the spacetime connection, 
are promoted to quantum operators. Then, the simplicity constraints are imposed at the quantum 
level as restrictions on group-theoretical data. The main challenge in this approach is essentially 
to impose consistently these second class simplicity constraints, and several proposals defining the 
different available spin foam models have been put forward in order to do so. 

The three spin foam models for four-dimensional quantum gravity that have been studied the 
most are the Barrett-Crane (BC) model [20], the Engle-Pereira-Rovelli-Livine (EPRL) model \21\ - 
123] . and the Freidel-Krasnov (FK) model |241 125j. They all rely on the general strategy outlined 
above, even though their construction requires to impose the simplicity constraints in drastically 
different ways. In the BC model, which was the first four-dimensional model to be introduced, the 
simplicity constraints are imposed as strong operator relations. This has the result of assigning 
simple representations of the gauge group to faces of the dual triangulation A*, and a specific 
unique intertwiner, known as the BC intertwiner |20l I26j . to edges of A*. Due to difficulties in 
reproducing the correct semiclassical limit [27H29] and the structure of the graviton propagator [30] . 
this model was partly discarded^ , and the search for new models was motivated with the additional 
hope to relate the spin foam quantization to the canonical structure of loop quantum gravity. This 
search culminated with the introduction of the EPRL and FK models, which are both based on 
the introduction of a linear version of the (originally quadratic) simplicity constraints [2^4j, and 
take as an important additional input the inclusion of the Barbero-Immirzi parameter. In these 
two models, the over-imposition of the simplicity constraints which is thought to be responsible for 
the issues with the BC vertex, is cured by using a weak imposition by means of coherent states for 
the FK model, or the Gupta-Bleuler scheme for the EPRL model. Our goal here is not to review 
extensively the details of these constructions, but rather to focus on the implementation of the 
simplicity constraints by introducing a model which allows for an immediate comparison between 
the alternative schemes. In fact, it has been argued on several occasions |33H35j that the current 
spin foam models miss the imposition of additional secondary second class constraints that are 
generated by the usual simplicity constraints. The usual point of view in the derivation of spin 
foam models is that it might be sufficient to impose the primary simplicity constraints consistently 
at all times, since the secondary constraints arise in fact as a consequence of this requirement. But 
even if this is true at the classical level, at the quantum level the secondary constraints are allowed 
to have non-vanishing fluctuations, and ought therefore be imposed as well. It is however very hard 
in the four-dimensional case to implement concretely this idea. In other words, we understand that 
there is a difficulty which needs to be resolved in the current four-dimensional spin foam models, 
but even if we have generic arguments indicating how to do so, a concrete realization has never 
been put forward. 

Let us summarize the situation. At the classical level, it is clear that the Plebanski action is 
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equivalent to the Palatini action once the B field is forced to be simple. But since we do not know 
how to compute the simplicial path integral for the Palatini action, we try to compute that of the 
Plebanski action by first quantizing the unconstrained topological theory and then imposing the 
simplicity constraints at the quantum level. In fact, it is quite clear that we are restricted in our 
understanding of how to properly perform this last step, by the fact that we do no know what to 
"expect" from the spin foam quantization of the four-dimensional Plebanski action. For this reason, 
it would be nice to have a model in which we could work out explicitly the two alternative methods. 
By this, we mean that we would like to perform the spin foam quantization of the action in which 
the simplicity constraints have already been solved at the classical level, and then compare this 
result to the spin foam quantization in which the simplicity constraints are solved in the quantum 
theory. This would be a way to test the various proposals to deal with the simplicity constraints. 
It has already been argued in the literature [36] that dealing with second class constraints at the 
quantum level can lead to inconsistencies, and is not compatible with the quantization program a 
la Dirac. However, the argument alone is not fully convincing, since it was formulated on a finite 
dimensional model with no clear analogy with gravity, and since it does not give a clear explanation 
for how the various impositions of the constraints in the four-dimensional spin foam models should 
lead to a model with inconsistent physical predictions. It might very well be that despite the fact 
that we are dealing with second class constraints at the quantum level, there exists a preferred 
scheme in which the simplicity constraints can be imposed in a (yet-to-be defined) robust way to 
lead to an acceptable model for quantum gravity. In the end, the only way to discriminate between 
the various proposals is either to extract physical predictions to compare with experiments, or 
at least to test the strategies on toy models which bear a close analogies with gravity. Let us 
also point out that a model was introduced in [33] in order to illustrate the claim that additional 
secondary second class constraints should be taken into account in the spin foam models. This 
model is based on the idea of reducing a four-dimensional S0(4) BF theory to an SU(2) one by 
means of simplicity constraints, and demonstrates that this is only possible if certain secondary 
second class constraints are imposed. The weak point of this model is that the constraints that it 
imposes are not derived from the Hamiltonian analysis of a given Plebanski theory, but are instead 
put in by hand. 

The aim of this paper is to formulate a robust model which allows to test the imposition of 
the spin foam simplicity constraints. Since pure gravity is always topological in three spacetime 
dimensions, it makes a priori no sense to have a Plebanski formulation in which the B field is 
forced to be simple. In fact three-dimensional BF theory is already equivalent to gravity, when B 
is a one- form field valued in the Lie algebra of SU(2). However, if one replaces the gauge group 
SU(2) by S0(4) (or its double cover Spin(4)), the three-dimensional theory admits simplicity 
constraints C similar to that of the four-dimensional theory. In the gravitational sector of solution 
to these simplicity constraint, the action becomes that of first order Riemannian gravity. In the 
topological sector, the theory is trivial. Therefore, at the classical level, the Plebanski theory that 
we introduce reduces to the Hilbert-Palatini action, whose spin foam quantization naturally leads 
to the Ponzano-Regge amplitudes. Now, it is also possible to perform the spin foam quantization of 
the three-dimensional SO (4) BF theory, and to impose the simplicity constraints at the quantum 
level, mimicking in particular the prescription of the BC and EPRL^ models. This allows for a 
direct verification of whether any of these proposals leads to the proper vertex amplitude. Our 
construction can be summarized in the diagram of figure [TJ 

We start in section |l] by introducing the classical framework for our analysis. In particular. 



^ Strictly speaking, we are going to study the EPR model [21] because we do not consider for the moment the 
inclusion of the Barbero-Immirzi parameter. 
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FIG. 1: The starting point of this investigation is the three-dimensional action Sbf- After imposition of the 
simpUcity constraints C, it reduces to the usual Hilbert-Palatini action, whose discretized partition function 
-Z^grav can be computed explicitly. The question we want to ask is whether it is possible to derive the same 
vertex amplitude by imposing the simplicity constraints at the quantum level on the partition function 2bf 
for the unconstrained action. In other words, does the imposition of the simplicity constraints a la BC or 
EPRL render this diagram commutative? If not, what is the additional input which is necessary in order to 
reconstruct the Ponzano-Regge model from the spin foam quantization of the three-dimensional Plebanski 
action? 



we show how the three-dimensional S0(4) BF action is reduced by the simplicity constraints to 
the usual action for three-dimensional gravity, and we perform the Hamiltonian analysis of the 
Plebanski theory to support this fact. Section |ll] is devoted to the study of the quantum theory. 
We first recall the state sum model corresponding to the constrained classical action, and then 
perform the spin form quantization of the Plebanski theory. For this, we start by writing down the 
simplicial path integral for the S0(4) BF theory, and then study the imposition of the simplicity 
constraints. In particular, we study the BC and EPRL prescriptions, in which the constraints 
are imposed strongly and weakly (in a sense that we will make more precise), respectively. We 
also study the so-called Warsaw modification |37ft39j of the EPRL model in order to illustrate its 
signification on the three-dimensional model at hand. Our conclusion is that the BC prescription 
cannot lead to the proper vertex amplitude, whereas the EPRL prescription, when supplemented 
with the imposition of the additional secondary second class constraints, enables one to recover 
the expected result. This might be a strong indication of the fact that a weak imposition of 
the simplicity constraints in the quantum theory is an important ingredient in order to derive the 
vertex amplitude of four-dimensional quantum gravity. Nevertheless, we want to emphasize that in 
our three-dimensional model, none of above-mentioned proposals lead by themselves to the proper 
vertex amplitude, which is that of the Ponzano-Regge model. The reason for this is clear: the 
current spin foam models miss the imposition of the secondary second class constraints that are 
exhibited in the canonical analysis of the Plebanski theory. Indeed, we show how these constraints 
can be solved by imposing conditions on the spatial components of the connection, and how their 
imposition in the path integral with the EPRL prescription allows to recover the proper vertex 
amplitude. 



I. CLASSICAL THEORY ACTIONS FOR THREE-DIMENSIONAL GRAVITY 

To introduce our model, let us start with the familiar BF action, written in three dimensions and 
for the gauge group S0(4). In terms of an so(4)-valued one-form field B^ (our notations are 
defined in appendix |a]) and the spacetime connection coj/ , it is given by 

Sbf[B,lo] = 1 [ d\e^'''PTr{B^,F,,) = \ [ d^x s^^-^pB^uFI^ , (1.1) 
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where Tr(- , •) denotes a trace in the Lie algebra so (4), and F^y is the curvature two- form associated 
to the connection. Clearly, this action is not that of Riemannian three-dimensional first order 
gravity, but this latter (or at least two copies of it) can be obtained if we use a decomposition of 
the fields into self-dual and anti self-dual components. Here we are not interested in doing this 
decomposition, because our goal is to mimic the constrained four-dimensional action which is the 
starting point for the spin foam quantization. To see that this is indeed possible, let us proceed in 
two steps. First, let us assume that the B field is of the form 

klX e^, 

where is a zero- form in M^, and a Lie algebra- valued one- form. This is the three-dimensional 
analogue of the notion of simple two-form field that is used in four-dimensional Plebanski theory. 
Second, let us choose a gauge (which we will refer to as the time-gauge) in which = (x") 0, 0, 0). 
This gauge choice reduces the symmetry group S0(4) to its SU(2) subgroup. A key observation is 
that it is always possible to fix x'^ = 1 due to the presence of the "scaling" symmetry 

x' ax', -ei, (L2) 



for a 7^ 0. With this choice of gauge, the action (1.1) becomes the usual action for gravity' 



SBF[B{e),u;] = %av[e,c^] = \ f d^x e'^^^^e^.F^^ (1.3) 

We have included the first equality to make explicit the fact that the three-dimensional Hilbert- 
Palatini action is simply a constrained BF action. This argument suggests that it is possible to 



obtain the action for gravity starting from (1.1) if we impose a simplicity- like condition on the B 
field, and pick up a rotational SU(2) subgroup of S0(4) by making a gauge choice for x^ (notice 
that this choice is not canonical). 

We are going to see in this section that the simplicity of the B field can be enforced with a 
Plebanski action whose Hamiltonian analysis will be detailed. 



A. Plebanski action and simplicity constraints 

As we mentioned above, the first things that we would like to do is to impose the simplicity of 
the B field in the S0(4) BF theory. Following what is done in four dimensions, we are going to 



introduce the Plebanski action as the sum of the BF action (1.1) plus a set of constraints on the B 
field. Our aim is to reproduce simplicity-like conditions ensuring that we are indeed dealing with 
three-dimensional gravity. The action that we consider is^ 

Spi[B, oj, <j)] = Sbf[B, oj] + constraints = J / d^x (e^'^f Tr(S^, F^^) + Tr(B^, *5^)) , (1.4) 



5 



This is in fact not straightforward, since the curvature in ( 1.3 1 involves connection components of the type ojji'. 
However, one can see that these terms appear quadratically as Lagrange multiphers, and that they are therefore 
vanishing by virtue of the equations of motion. 

Notice that the constraint term is not a three-form since the spacetime indices of the B field are contracted with 
those of the multiplier 0, and not those of the antisymmetric Levi-Civita symbol. This is similar to what happens 
in the four-dimensional Plebanski theory when the multiplier acts on the spacetime indices of the B field. 
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where ★ denotes the Hodge dual operator in so (4), and (j)^'^ is a symmetric Lagrange multipHer 
used to enforce the constraints 

C^, = Tr(5^, *B,) = ^eijKLBj/B^^ « 0. (1.5) 

We want to understand the exact meaning of these constraints in the present context. The key 
difference with four-dimensional Plebanski theory is that we are dealing here with vectors, and not 
bivectors. Therefore, we need another notion of simplicity than the usual one used in spin foam 
models. As we have already mentioned earlier, we are going to say that a vector Bj^^ is simple if 
it can be written as 

r>IJ _ A/^fl-^ 

klX e^, 



where is a zero- form in M^, and a one- form. The fact that the constraint C defined by (1.5) 
can be used to obtain simple vectors is ensured by the following proposition: 

Proposition 1. When the simplicity constraints C are satisfied, there are three possible solutions 
for the B field. They correspond respectively to: 

• A gravitational sector, in which there exists an M.'^ -valued one form and a vector in 
such that 



Bl' = e'^KLx'^ej:, Ki ^ B^; = (x A e,)\ ^ ^ ^.B^' = x% - x'e'. 



• A topological sector, in which there exists an M.^ -valued one form and a vector x^ in 
such that 

B'/ = ^e'^KLX^'e]: = x'ej, - x' Kj, = x'ej, - x'e% L], = (x A e^)\ 

• A degenerate sector, in which det(i^^) = and det(Lp = 0. 

Proof. Let us introduce the boost and rotational components of the field B^^^ . They are given 
respectively by K'^^ and L^. Each of these one- forms can be interpreted as a vector in M^, and the 
simplicity constraints simply mean that the vectors are orthogonal to the vectors Ly, i.e. 

L^-Ky = ^, V/x,z/. (1.6) 

Let us first assume that the three vectors form a basis of M^. Due to the simplicity constraints 
(1.6), the three remaining vectors K^^ are linked, and therefore they lie in the same plane, whose 
(non-necessarily unit) normal is denoted by x- More precisely, there exist three vectors such 
that 

i^M = X A e^. (1.7) 

Furthermore, since the vectors and K^^ are orthogonal, we can write to 

= x'e^ - elx, (1-8) 

where x^ is a (non-vanishing) scalar, and is a real-valued one form. Due to the fact that the 
vectors form a basis, the three vectors form a basis as well. The solution given by (1.7) and 
(1.8) corresponds to the gravitational sector. 
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If we assume on the contrary that the three vectors form a basis, the same construction can 
be used to obtain the topological sector described in proposition [TJ 

When det(-fC^) = and det(L^) = 0, i.e. neither the nor the form a basis, we say that 
the B field belongs to the degenerate sector. ■ 

It is clear from the expression of Bj/ in the gravitational sector that in the time gauge, where we 
have X* = and x'^ = 1, the boost component Kj^ is vanishing, while the rotational component 
reduces to the triad e^. This supports the heuristic argument that we have given at the beginning 



of this section in order to derive the action (1.3) from the SO (4) BF theory. Note also that in the 
gravitational sector we have det(i^^) = 0, while in the topological sector we have det(L^) = 0. 

Let us conclude with a remark. There are 6 simplicity constraints C^jy acting on the 18 com- 
ponents Bj^^ of the B field. Therefore we expect the simple B field to be written in terms of 12 
components only. Here, we have expressed B in terms of the 16 components ^-^id e^. However, 
notice that there are 4 redundant components due to the presence of the following symmetries: 



A rescaling symmetry, given by (1.2), which allows to remove one component (for instance, 
we can fix x^Xi = !)• 

Three translational symmetries acting in the non- degenerate sector. They are generated by 
a vector and act like 



As a consequence, we can always choose = 0, but here we do not make this choice in 
order not to break the SO (4) symmetry of the theory. 

B. Gravitational and topological sectors 

Now that we have a good understanding of how the simplicity constraints can be used to define 
a sector of the Plebanski theory in which the vector Bj/ is simple, let us see how it is possible 
to recover the action for three-dimensional gravity. For this, we are going to study separately the 
gravitational and topological sectors introduced in proposition ([T]). 

Using the decomposition of so (4) into self-dual and anti self-dual generators (see appendix [B|) , 
we can write the B field as 



where ^"""'-Bi = Tr(i3^, *='='J*), and the BF action (1.1) therefore becomes 



SbfI^^^B, ^-^B, = \j^ d3xe^''^(Tr((+)i3,„ ^+^F,p) + Tr((-)i?^, ("'F.p)) . (1.9) 

As usual, the equations of motion with respect to the connection ''^^oj lead to the torsion-free 
condition 

r((±'5, = 0, 

and if dei{^'^^B) ^ 0, this relation can be inverted to find the torsion-free spin connection. This 



latter, when plugged back into the original action (1.9), leads to the sum of two second order 
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Einstein-Hilbert actions, 

^EH^V/'V] = / d'x^/FV\n^^'9^u) + le- [ d3x^^7^(''V)> (1-10) 
^ Jm ^ JM 

each being defined with respect to a two-dimensional Urbantke-like metric 00] ^^''g^u = ^^^Bfj_ ■ 

(in the sense that it is constructed with the B field). In this expression, denotes the sign of 

det('*'i?) (we refer the reader to appendix [c] for the calculation of e^). 

In the gravitational and topological sectors, the self-dual and anti self-dual components of the 
B field are given by: 

• gravitational sector: = =F(x A e^)* -|- (x*e° — x°ep, £"•" = £", (1.11a) 

• topological sector: = — (x A e^)* it (x*e° — x^e^), e~^ = —e~. (1.11b) 

In each of these two sectors, we can now compute the two Urbantke metrics '-^^g^.u defined by the 
expressions above. In fact, from the simplicity constraints, we know already that ^^''g^v = '"~^g^u- 
It is however less trivial to see that the metrics in the topological and gravitational sector are 
identical. A simple calculation shows that we have 

= V = 9t.u = (e^ • e.)(x' + (X°)') - (x • e^)(x • e.) + x'e°e° - x\ ■ (e°e. + e^e^), (1.12) 

where x^ = X^Xi- Iii the time gauge, this metric reduces to g^u = e-u- 

Gathering these results on the Urbantke metric and the sign factors e^, we can conclude that 



the action (1.1) reduces in the gravitational sector to 

SFRig^iu] = / d^x^J\g\l^(g^y), 
JM 

while in the topological sector it becomes simply S-E\^\g^i\ = 0. 

This shows that the gravitational sector to the solutions of the simplicity constraints C corre- 
sponds indeed to three-dimensional gravity, whereas in the topological sector the action vanishes 
on-shell. This is exactly what happens in the four-dimensional Plebanski theory. Interestingly, the 
presence of a topological sector in this three-dimensional theory allows for the introduction of a 



Barbero-Immirzi parameter, but we shall come back to this point in section III 



C. Canonical analysis of the Plebanski action 

In this subsection we perform the canonical analysis of the three-dimensional SO (4) Plebanski 
theory. It is similar to the study of the four-dimensional Plebanski action [41]. Using a 2 -|- 1 



decomposition of the spacetime manifold, (1.4) becomes 

5pi[S,u;,(/.] = j^dtj^d^x (^-e^'Ti{Ba,douib)+le'''i:tiBo,Fab)+e'''Tricjo,Tab)+lr'^^ 

where the curvature and torsion two-forms are defined by^ 

e'^'Fab = e''\daUJb - dbCOa + [iOa,U:b]), e^'Tab = e^\daBb + [Ua, Bb]). (1.13) 



^ The torsion is here defined up to a factor 1/2 for convenience. Its exact expression is Tab ~ daBt — dbBa + [lOa, Bb\ 



[U)b, Ba]- 
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Wc can see from this formula for the action that ujq, Bq and (p^'^ are non-dynamical variables, since 
the Lagrangian does not feature their time derivatives. The time component coq of the connec- 
tion and (f)^'^ are true Lagrange multipliers, enforcing respectively the torsion-free and simplicity 
constraints. However, since Bq is involved quadratically in the simplicity constraint, it cannot be 
treated as a Lagrange multiplier. For this reason, we add to the Lagrangian the term 

Tr(7ro, doBo) + Tr(^o, ttq), 

where ttq and fiQ are new auxiliary S0(4)-valued fields that do not affect the dynamics of the theory. 

The basic variables of the theory are then the 12 spatial components w^"^ of the connection, 
their 12 canonical momenta -B^"^, and the 6 components Bq"^ with their momenta ttq"^ . The Poisson 
structure is given by 

{5,"(x),a;f ^(y)} = -S'^'^^Sab6\x - y), {i?o"(x), 7ro^^(y)} = d''''''^6\x - y), (1.14) 
where 6^-^'^^ = {6^^ 6-^^ - 6^^6'^^)/2. The total Hamiltonian is 

H = -J^d^x Qe"^ Tr(So, F^b) + e"' Tr(a;o, Tab) + Tr(i?^, *B,) + Tr(//o, ttq) j . (1.15) 

We can now identify the primary constraints and compute their Poisson bracket with this Hamil- 
tonian in order to study their evolution in time. 

1. Primary and secondary constraints 
We have the following 18 primary constraints: 

Cf,^ = TV(5^, *B^) « 0, Tab ~ 0, TTo « 0. 

The first set is obtained by varying the action with respect to <p'^'^, the second set by varying the 
action with respect to ljq, and the third one by varying with respect to /ig 

Before going any further, let us introduce smeared variables in order to deal with the delta 
functions appearing in the Poisson brackets. To this end, we define the quantities 

T{u)= [ d'xe''^Tv{u,Tab), F{u) = [ d'x e"'' Tj:{u, Fab), 

for any S0(4)-valued smooth test function u{x) on S. This enables us to compute 

{F{u),F{v)}=0, {Tiu),T{v)} = -T{[u,v]), {T{u), F{v)} = -Fi[u,v]), 

and 

{T{u),Ba} = [u,Ba], {T{u),COa} = -T^aU, {F{u), Ba} = -2VaU, {F{u) , OJa} = 0, 

where T> stands for the covariant derivative with respect to the connection w. Now we can compute 
the Poisson bracket of the primary constraints with the total Hamiltonian, to see if the requirement 
that they be preserved under the time evolution gives rise to secondary constraints. 
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a. Evolution of the constraints Tab 



First, it is easy to show that the conservation of Tab does not lead to any secondary constraints. 
To see this, we can introduce the new constraint 

G(n) = T{u) - d^xTiiu, [B^M), (1-16) 

and notice that it is in fact the generator of the infinitesimal SO (4) gauge transformations on the 
phase space. Indeed, on the action of G{u) on the phase space variables (1.14) is given by 

{G{u),B^] = [u,B^], {G{u),UJa] = -Vau. 

Since the Hamiltonian contains only terms with traces, it is left invarian t und er the action of this 
constraint, which therefore satisfies G = {if, G} ~ 0. The modification ( |l.l6[ ) of the constraint is 
permissible since it amounts to adding a term which is itself proportional to a constraint. 



h. Evolution of the constraints ttq 



The requirement that the primary constraint vro be preserved in time leads to the following relation: 



Among these 6 equations, 3 are in fact fixing the components cjf'^ (if we assume that the non- 
degeneracy condition holds), and the remaining 3 have to be imposed as secondary constraints. 
Projecting onto the vector Z?^, these secondary constraints can be written as 

*^ = e"''Tr(F,5,i?^)~0. 



c. Evolution of the simplicity constraints C^, 



To study the evolution of the simplicity constraint, it is convenient to compute separately the 
Poisson bracket of its various components with the Hamiltonian. Requiring that the constraints 
Coo and Coa be preserved under time evolution leads to 

Coo = -2Tr(/xo,*5o)«0, (1.17a) 
Coa = Tv{*Bo,VaBo) + Tr(*So, [Ba,UJo]) - Tv{fio,*Ba) « 0. (1.17b) 

These 3 equations determine 3 of the 6 components of the multiplier ^O; and imply no secondary 
constraints on the dynamical variables of the phase space. Finally, a direct computation shows 
that the requirement that Cab ~ leads to the 3 secondary constraints 

= TiiVaBo,*Bb) + TT{VbBo,^Ba) « 0. (1.18) 



The Dirac algorithm stops here, and there are no tertiary constraints. 
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d. Summary of the analysis of the primary constraints 

To summarize, we have a theory with 18 primary constraints consisting of 6 constraints ttq, 6 
constraints Tab, and 6 constraints C^u- They generate the 6 secondary constraints comprising the 
3 constraints and the 3 constraints 



2. First and second class constraints 

In the previous subsection, we have derived the primary and secondary constraints of the Plebanski 
theory. We are now going to spht them between first class and second class constraints. 



a. The first class constraints 



First, it is easy to see that the constraints G are first class. Indeed, they commute with all the 
other constraints since they generate the infinitesimal S0(4) gauge symmetries. 

The analysis of equation (1.17) suggests that amongst the 6 constraints vro ~ 0, 3 are first class 



and 3 are second class. To see that this is indeed the case, let us decompose the set of constraints 
ttq ~ into 



fC^, = Tr(7ro, 



/C^ = Tr(7ro,*-B^). 



If the B field does not belong to the degenerate sector, the previous set of constraints is equivalent 
to the requirement that ttq — 0. Furthermore, a direct computation (see the algebra of constraints 
in appendix IdI) shows that the constraints JC^ are first class. They are associated with the 3 



components of fiQ that remain undetermined after taking (1.17) into account. 

Finally, the last first class constraints are given by a new set which is obtained by adding to 
a linear combination of constraints. In particular, is defined as being equal (up to a factor 



2) to the Hamiltonian density (1.15): 

§0 = ^0 + 2e'^''Tr(wo,Tafe) + cj)'"' Tr{Bf„*B,) + 2Tr(^o,vro), 

where the values of the Lagrange multipliers ojq, (p^'^ and /xq, are those determined by the constraint 
analysis. When the Lagrange multipliers are not fixed by the Hamiltonian analysis, it is possible 
to fix them to the value zero in the expression for ^'q. Concerning the constraints ^'a, a direct 
calculation shows that they are given by 

$a = + 2e^' Tl{uJa, Uc) " Tr(^o, daBo) 

= 2e'" {Tr{Ba,dbOOc) + TV(w„ ObB,)) - Tr(7ro, daBo), 

and that they generate as expected space diffeomorphisms on the phase space variables: 

{§(e),Wa} = -Ci:COa, {5(0,^/.} = 

where ^ is a vector field on S, and ^(^) denotes the smearing of the ^a with 



12 



b. The second class constraints 

The remaining constraints, 0^,^, and K^, are the second class constraints of the theory. In 
order to prove this, it is possible to compute the (square 12-dimensional) Dirac matrix A, whose 
elements are given by Poisson brackets between the various (candidate) second class constraints, 
and show that its determinant is non-vanishing. This is indeed the case, since the Dirac matrix is 
given by 







Cab 




^ab 


Coi/ 


= 


= 




= 


Ccd 


= 


= 


= 


Mab,cd 




-Km 


= 


« 


Nab,u 




= 


—Mcd,ab 


-Ncd,t. 


Pab,cd 



and it satisfies clearly det(A) ^ 0. Here the elements of the matrix V are given by 
V = {^M,Co.} = {'IV(7ro,*5^),IV(So,*B,)} = W^^, + = 2^^,, 

and we have |det(y)| oc V^, where V = ^ det((?^;y) denotes the three-dimensional volume. The 
matrix M has elements determined by the Poisson brackets 

Mab,cd = {$a6,Ce4 = ecaTr([5o, ^d, ^d) + £daTr([5o, B^lB,) 

+eefeTV([5o, Ba\,Bd) + edbTV([5o, ^a], ^c), 

and satisfies |det(M)| oc V^. The explicit form of the matrices Nab,v = {"^afe^^j^} a-nd Pab,cb = 
{$aft, ^cd} is not important in order to prove the invertibility of the Dirac matrix, and it is in fact 
easy to see that its determinant is given by^ 

det(A) = (det(l^)det(M))^ = cV^°. 
3. Summary 

To summarize, the first class constraints of the system are given by /C^, ^'^^ and Gab, £^nd the 
constraints $a6 and K^, are of second class. Therefore, we have 36 phase space variables B^ 
and Loj^^ that are subject to 12 first class constraints (generating gauge symmetries) and 12 second 
class constraints, and the usual counting shows that there are no phase space degrees of freedom. 
This is of course to be expected in three-dimensional gravity. 

As usual, the first class constraints of the theory are the infinitesimal generators of the gauge 
symmetries. The constraint G is the generator of the internal S0(4) symmetries. The 3 primary 
constraints /C^^ appear because we have treated the lapse and the shift encoded in B as dynamical 
variables. The constraints /Co and ICa therefore encode the vanishing of the momenta ttat and ttat" 
conjugated to the lapse and the shift. The remaining constraints, *o and ^'a, are related to the 
scalar and vector constraints, i.e. to the spacctime diffeomorphisms. 

Let us conclude this canonical analysis with an important remark concerning the second class 
constraint. The 3 constraints K/j, constrain 3 out of the 6 components of the variables ttq"^ , which 
are the conjugate momenta to B^'^ . As a consequence, only 3 components of tto"^ are left, and they 



The exact proportionality coefficient is c = 2^^3®, since | det(F)| = 2^V^, and | det(M)| = 2"3^V^. 
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correspond to the momenta vrjv and vrjva associated with the lapse and shift variables. Moreover, 
the lapse and the shift are the 3 independent components of B^^ that are left after imposing the 
3 components Cq^ of the simplicity constraints. The spatial part Cab of the simplicity constraints 
ensures that the components B^^ can be expressed in terms of the zero-form and the one-form 
e^. Finally, the meaning of the 3 constraints ^^ab is that 3 of the components of the connection 
u:l^ can be expressed in terms of and ■ It is important to stress that these are secondary 
constraints which have been obtained from the requirement that Cab be preserved in time. In 
other words, the usual simplicity constraints are of second class because there are in fact secondary 
constraints which do not commute with the primary ones. This point will be very important in the 
spin foam quantization that we perform in the next section. We are going to show the importance 
of imposing the secondary second class constraints in the spin foam models. 



II. SPIN FOAM QUANTIZATION 

In this section, we study the quantization of the three-dimensional Plebanski theory. From the 



classical analysis of the previous section, it should be clear that the action (1.4) is equivalent in 
the gravitational sector to that of first order gravity, whose spin foam quantization is the Ponzano- 
Regge model. We can now try to reproduce the strategy of four-dimensional spin foam models, and 
see if the spin foam quantization of the full Plebanski theory leads to the Ponzano-Regge model 
as well. This is a well-posed question which will enable us to clarify the issue of imposing the 
simplicity constraints at the quantum level. 



A. The Ponzano-Regge model 

We have seen in the previous section that when the simplicity constraints are imposed at the 



classical level, the gravitational sector of solutions to the Plebanski action (1.4) is equivalent to 
the SU(2) Palatini action. It is therefore possible to compute explicitly the simplicial path integral 
of this theory, and we know that it leads to the Ponzano-Regge state sum model, whose vertex 
amplitude is given by an SU(2) Qj symbol. We recall some basic facts about the derivation of this 
result, and also introduce the notations that we will use in the next subsections. 

The partition function of three-dimensional Riemannian gravity is formally given by 



^grav = j d[e\d[uj] exp [i d^x e^'^P tr(e^, F^p)) 



where tr(- , •) denotes a trace in the Lie algebra su(2). To compute the simplicial path integral, 
we introduce a simplicial decomposition A of the spacetime manifold along with its dual two- 
complex A* |42j . This latter is consisting of vertices v (dual to tetrahedra r G A), edges e (dual 
to triangles i G A) and faces / (dual to links £ G A). Since the triad ej^ is a one- form, it is natural 
to integrate it along the one-cells (links) of A to obtain su(2)-valued elements Xf. The connection 
is discretized by computing its holonomy he along the edges e of A*. We can then discretize the 
curvature by taking the product of holonomies along the edges lying on the boundary of a face 
/ G A* to form 



^/ = n 

ecf 
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With these variables, the discretized version of the action (1.3) can be written as 



SBF[Xf,hf]= tr{Xf,hf), 
/6A* 



and the partition function becomes 



^grav = ( n / I n / d/le I exp i V tT{Xf, hf)], (2.1) 

where dXf is the Lebesgue measure on su(2) ~ M'^, and d/ig the Haar measure on SU(2). It is now 
possible to perform the integral over ^/ to obtain 



/ d/le TT (^SU(2)(^/), 

. JSU{2) j j.^^. 



where the delta distribution over SU(2) imposes the flatness of the connection. Using the Peter- 
Weyl decomposition, we can write 

= E n (2^/ + 1) ( n ) n x., ( n ) , (2.2) 

{i}^{/}/eA- yeeA*-'^U(2) J feA* \eCf j 

where the sum is taken over all the possible SU(2) representations j labeling the set of faces / G A*. 
For arbitrary cellular decompositions A, let us call n the number of faces meeting at every edge 



e E A*. In (2.2) we will have an integral over hf, of n products of representation matrices for the 



group element /ig- We can therefore use the fact that 

/ d/leD(^i)(/le)...D(j")(/ie) = ^iei: (2.3) 

projects onto the space Invgu(2) (^su(2) ® " ' ® ^su(2)) °^ intertwiners between the n representa- 
tions coloring the faces meeting at the edge e. Then, all the intertwiners meeting at a vertex v can 
be contracted to define a vertex amplitude ^^((jj-dd, ^cdu)- Finally, the partition function can be 
written as a sum over spin foam amplitudes 

^grav = ^ ^ Wi'^jf + ^) W A^{3f^v,^e^v)■ 
To clarify the meaning of this formula, let us assume that the cellular decomposition A is simplicial. 
In this case, vertices in A* are four-valent, while edges are three-valent. The vertex amplitude is 
therefore given by a contraction of four three-valent intertwiners, which is the so-called 6j symbol. 
The final result of this calculation is 

^grav = E n + 1) n {6J} (2.4) 

{i}^{/} /eA* 1-eA* 
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where the labels ji, . . . j'e are associated to the six faces dual to the links I a tetrahedron r G A. 
Notice that the sum over intertwiners has now disappeared since there is only a unique (up to 
normalization) three-valent intertwiner. In other words, the intertwiner space 



Invq 



SU(2) I ^SU(2) "•SU(2) 



(ii) 



ih) 



SU(2) 



is one- dimensional when (ji,j2)j3) satisfy the triangular inequalities, and zero-dimensional other- 



wise. 



It is convenient at this point to introduce the diagrammatic notations [l3] , in which lines are used 
to represent unitary irreducible representations, and boxes to represent integrals over the group 



defining a projector on the intertwiner space following (2.3). With this notation, the partition 



function defining the Ponzano-Regge model can be represented as follows: 



I I I 




TTT 



Once the integrations are performed, we obtain a contraction a four three-valent intertwiners for 



each vertex dual to a tetrahedron, and the partition function becomes (2.4) 



B. Quantization of the SO (4) BF sector 

Spin foam models in four dimensions take as a starting point the four-dimensional Plebanski 
action, and are derived by first quantizing the topological part of the action, and then imposing 
the simplicity constraints at the quantum level. 

Here we follow the same procedure, and start by quantizing the S0(4) BF part of the Plebanski 



action ( 1.4 ). This can be done very easily along the lines of the standard construction introduced in 
the previous subsection, with the only difference that SU(2) is now replaced by S0(4). Therefore, 
the B field is discretized by assigning an 50 (4) element Bf to each face / S A* dual to a one-cell 
in A. Once we perform the integration over Bf in the simplicial path integral, we are left with 

2bF= n / d/le TT (5so(4)(/i/), 

where now we have a product of delta distributions on the group S0(4), which can again be 
evaluated with the Peter- Weyl decomposition. Using the fact that S0(4) ~ SU(2) x SU(2)/Z2, it 
is possible to write an irreducible representation p of S0(4) as a tensor product p = j+ of two 
irreducible representations of SU(2). The discretized partition function for the three-dimensional 
SO (4) BF theory then takes the form 

^BF = E n (2^; + 1)(2j7 + 1) n {6i+}{6j-}. 
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Now, we have to implement the simplicity constraints in this state sum in order to ensure that it 
describes quantum three dimensional gravity. Indeed, since at the classical level the gravitational 
sector of the Plebanski action corresponds to the action for gravity, whose spin foam quantization 
corresponds to (2.4), we expect that a proper imposition of the simplicity constraints will lead to 
the same result. 



C. Imposition of the simplicity constraints 

In the simplicial path integral that we have written for the three-dimensional SO (4) BF theory, 
the one- forms Bj/ have been discretized by assigning an element By E 50 (4) to the faces of A*, 
which are dual to the links i G A defining the boundary dt of the two-cells (triangles) t € A. These 
vectors in fact determine the geometry of the triangles, and they satisfy the closure constraint 

E = 0- 

f&dt 

This is just the discrete analogue of the continuous Gauss law which ensures gauge invariance under 
the action of SO (4). Additionally, the discretized B field is required to satisfy the discretized version 
of the simplicity constraints C, which are given by^ 

diagonal simplicity: euKiByBf^ ^0, ^ f e dt, (2.5) 

off-diagonal simplicity: eijklBY Bf^ k, yfj'Gdt. (2.6) 

These simplicity constraints do not distinguish between the topological and gravitational sectors, 
since they are left unchanged if we change B for -kB. 

Just like in the four-dimensional theory |24| I35j. the simplicity constraints can be linearized, 
and are equivalent to the requirement that 

xiBy = o, yfedt. (2.7) 

Indeed, in the gravitational sector where we have B^ = e^'/^/^X^^^, it is clear that the linear 
simplicity constraint is satisfied. This linearized version has therefore the advantage of selecting 
the gravitational sector. 

The quantization of the theory is based on the symplectic structure of the topological part of 
the action. In the BF action, the variable B is canonically conjugated to the connection w, and 
the quantization rule is therefore simply to identify the discretized field Bf with the generators 
13 f = Jf of the Lie algebra so (4). The next step is to impose the simplicity constraints on these 
generators following the different spin foam models that have been introduced in the literature. 

The imposition of the constraints is done before the integration over the connection components 
defining the intertwiners. Schematically, we can represent this imposition of the constraints as a 
white box acting on the group generators between two neighboring triangles. Our task is therefore 



We do not mention the additional volume constraints, since they are automatically implemented when the closure, 
diagonal and off-diagonal simplicity constraints are satisfied. 
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to implement the simplicity constraints in such a way that the partition function reduces to 



{i+J-}^{/}/6A- 





1 




1 


L 1 










i 










1 — 



Here we have represented the implementation of the constraints via C on one triangle only for the 
sake of clarity. More explicitly the meaning of this graphical notation is that we have an integral 
over the S0(4) holonomies, together with a yet-to-be defined implementation of the simplicity 
constraints C. In other words, we have 

E = E = j dhS{C), (2.9) 

and the implementation of the simplicity constraints is done before the integration over the 
holonomies that defines the intertwiners leaving on the edges. 



1. The BC prescription 

We are going to start by solving the diagonal simplicity constraint. It is known that this constraint 
is equivalent to the requirement that the pseudo-scalar quadratic Casimir operator of so(4) defined 
by C2 = Tr(*J, J) vanishes, which constrains the representations to be simple. Indeed, in terms of 



the self-dual and anti self-dual generators, (2.5) takes the form 



Tr(*J/, Jf) = Tr((<+)J/ - <-'J/), C+'J/ + ("U/)) = Lf ■ Kf = 0, 

which implies the restriction = jj = jj to simple representations of S0(4) |20l I44j . The off- 
diagonal simplicity constraint involves two different faces / and /' dual to links i and i' belonging 
to the same triangle. They are constraints on the intertwiner states of the quantum theory, and if 
we impose them strongly we obtain the unique Barrett-Crane intertwiner (2()j . which selects only 
the subspace jj = of the decomposition 



n 



SO(4) 



2i/ 



0/) 

SU(2)- 



(2.10) 



Indeed, it is clear that the strong imposition of the constraint Ci = Tr{Jf,Jf) = 2jf{jf + 1) 
imposes jf = 0. Therefore, the operator implementing the simplicity constraint in the BC model 
acts like 



^Bc(jl, j2, js) 



(2.11) 
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As a result, if we use this operator in the path integral (2.8), we obtain the state sum 



1 1 1 1 1 1 



{3}^{f} 




I I I I I I 



E n (2^/ + n 



Notice that the BC operator defined in (2.11) is not a projector. However, it can be renormalized 
to give the genuine projector 



^jl ^33 



where dj. = 2ji + 1 denotes the dimension of the spin-j representation. This modification only 
affects the weight associated to the faces, and not the vertex amplitude. 

We see from this construction that the BC prescription, in which the simplicity constraints are 
imposed strongly, does not lead to the proper vertex amplitude. The next step is to check wether 
a weak imposition of the constraints can cure this discrepancy. 



2. The EPRL prescription 



Now we want to test the EPRL imposition of the simplicity constraints. The diagonal simplicity 
constraint can again be imposed strongly, and leads to the restriction = jJ = jf to simple 
representations. The idea is then to impose the off-diagonal constraint weakly. For this, we use 



the linear constraint (2.7), and notice that in the time gauge it becomes 



B 



K 



f 



0. 



One way to impose this constraint weakly is to use the master constraint trick [2T] , and impose the 
vanishing of the square of the boost components Xj. From the expression of the quadratic Casimir 

operators of S0(4), we find that this amounts to representing C by a projection and an integration 
over the highest weight representations, i.e. jJ + jJ = 2jf. The EPRL operator realizing this 
imposition of the simplicity constraints can be represented as 



^EPRL(jl, j2,j3) 



(2.12) 
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where the integration is taken over the spin-2jj representations. With this implementation of the 
simphcity constraints, the partition function (2.8) becomes 



1 1 1 1 1 1 



^EPRL = J] n + 1)' E 

{j}^{f} feA* 




= E n + n ( n f-^^-^n > 

with a, (3 = 1, . . . , 4, the index a labehng the edges, and the pair af3 labeling the faces. The fusion 
coefficient fa{ji, j2, js) attached to each edge e is given by the following 9j symbol: 




(2.13) 



Notice that the dependance of the fusion coefficient on the intertwiners i, i"*" and i~ , drops out 
because these intertwiners are three-valent and therefore unique. This is however an artifact 
of the fact that we have chosen the cellular decomposition A of the spacetime manifold to be 
simplicial (each edge of A* is therefore bounded by three faces). Another choice would have lead 
to intertwiners with higher valency. 



3. The Warsaw modification of the EPRL prescription 



Just like in the four-dimensional case, the operator (2.12) that we have in the previous subsection 



to implement the simplicity constraints is not a projector, i.e 

2 



E - 



E - 

P - 

R Z 

L - 



This "problem" can be solved following the idea introduced in the four-dimensional context in 
and defining a normalized operator 













E 






P 






R 






L 



E _ 

P - 

R = 

L - 



which is now a genuine projector. The 9j symbol appearing here is the one given by (2.13). This 
normalization factor is in some sense very trivial because we are dealing with a three-dimensional 
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theory and the sum over the (three- valent) intertwiners has disappeared. As a result, if we use 
this projector to impose the simphcity constraints, the 9j in its definition compensates the fusion 
coefficient obtained from the EPRL operator, and we obtain that the vertices of the BC and 
Warsaw-EPRL models are the same. 

Therefore, in the three-dimensional model at hand, and in the absence of Immirzi parameter, 
the prescription of defining a true projector to impose the simplicity constraints does not lead to 
the proper vertex amplitude either. 



D. Imposition of the secondary second class constraints 



In the previous subsection, we have tried to impose the simplicity constraints in the partition 
function for the SO (4) BF theory in order to recover the spin foam quantization of the gravitational 



sector of the classical Plebanski theory ( 1.4 ), i.e. the Ponzano-Regge model. The partition functions 



Z-Qc and Z-EPRL that we have obtained differ from that of the Ponzano-Regge model by the form 
of the face and vertex amplitudes. 

The form of the face amplitude is supposed to be related to the path integral measure. In the 



spin foam quantization (see for example equation (2.1)), the measures are taken to be the Haar 
measure for SU(2)- or S0(4)-valued elements, and the usual Lebesgue measure for Lie algebra- 
valued elements. However, it is known that the Plebanski action has a non-trivial path integral 
measure due to the presence of the various types of constraints [411 H5l H6] . It would be interesting 
to use the results derived during the canonical analysis of section |ll] (in particular, the determinant 
of the Dirac matrix) to compute the path integral measure of the Plebanski theory (1.4), and see 
if we can recover the face amplitude {2jf + 1) of the Ponzano-Regge model. We plan to come back 
to this question in future work. 

For the time being, let us focus on the vertex amplitude. To derive the spin foam models of the 
previous subsections, following the prescriptions of the BC and EPRL models, we have implemented 
the second class simplicity constraints C. However, as long emphasized by Alexandrov [HI [33H35] 
for instance, the primary simplicity constraints C generate secondary constraints with whom they 
form a second class pair, and these secondary second class constraints should also be imposed in 
the derivation of the spin foam models. Indeed, we have seen in the canonical analysis of the 
Plebanski action in section I C that the requirement that the primary simplicity constraint Cnu be 



preserved in time generates the secondary constraints ^ab- These constraints are additional second 
class constraints that have to be implemented somehow. 

Let us have a closer look at the constraints ^ab given by (1.18). For the sake of clarity, we 



only focus on the first term, which can be written with the decomposition into self-dual and anti 
self-dual fields in the form 

TriVaBo,^Bb) = TridaBo, *Bb) + Tr([c^„ Bq], *A) 

= Tr(9„(+'i?o, - Tr(a,(-'Bo, '-'Bb) 

+ Tr([(+W, '^'Bo], '+'i?fe) - Tr([(-W, '-'Bq], ("'S,). 



From equation (1.11a), we see that in the time gauge and in the gravitational sector we have 



^"•"^B^ = ' 'B^. The previous equation therefore reduces to 



Tr(P„Bo,*Bb) =Trl 



r(+), 



(-), 



and we have a similar relation for the symmetric term. It is therefore clear that the condition 
= '"'a; solves the constraint (in the gravitational sector and in the time gauge). To 
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implement this additional constraint at the level of the partition function, let us first discretize the 
condition ^+'a; = '"'a;, and write it in terms of holonomies as = We can then replace the 



operator (2.9) by 



dh6{<l>)6{C) = J dh5{^+^h^-^h-^)6{C), 



where the white box should implement the simplicity constraints C, and the integration measure 
has been changed to take into account the imposition of the secondary second class constraints ^. 

If we impose the simplicity constraints 6{C) a la BC, i.e. by projecting first onto simple repre- 
sentations = jj = jf, and then considering only the lowest weight module in the decomposition 



(2.10), we obtain a trivial vertex amplitude which is simply given by 



A — dj^ dj,^ dj^ dj^ dj^ . 

These weights can be seen as associated to the faces, and not to the vertices of A*. Therefore, 
there is no non-trivial dynamics associated with the vertices. 

Instead, if we use the weak imposition of the simplicity constraints, following the EPRL pre- 
scription, to compute and constrain the partition function of the SO (4) BF theory, we obtain the 
state sum 



{j}^{f} feA- 



{j}^{f} /6A* 



AAA 



T T T 



2ji 



2J3 
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AAA 



E npi/+i)=ns- 



T Y Y 

1 [ 2ji 2i2 2i3 

,76 I 2j4 2j5 2^6 



{i}^{/} /eA* 

The vertex amplitude is now given by a 6j symbol, which is what we expected. 



E. Discussion 



Let us make a few comments about the result obtained above. The first one concerns the weights 
associated with the faces. By taking into account the secondary second class constraints, we did 
not recover the correct weights for the faces. As we already said, this comes from the fact that we 
should have used the correct non-trivial path integral measure. It would be interesting to check 
wether taking into account this measure leads to the correct face amplitude. 

The second remark concerns the fact that the 6j symbol associated with the vertices involves 
the representations 2jf and not the representations jf. This might very well be a consequence of 
the fact that the state sum model that we have constructed corresponds classically to two copies 



of the Einstein-Hilbert action, according to (1.10) 



Finally, we would like to discuss the relevance of the EPRL prescription compared to the BC 
one. Both prescriptions are constraints on the intertwiner degrees of freedom labeling the edges 
of A*, while the secondary second class constraints that we have imposed act on the faces. As 
we have seen, even with a proper implementation of the secondary second class constraints, the 
BC prescription is too strong and "kills" all the dynamics. It is therefore clear that the simplicity 
constraints C should be imposed weakly, which is an argument in favor of the so-called new models. 



III. NOTE ON THE INCLUSION OF THE BARBERO-IMMIRZI PARAMETER 

We would like to end this work with a note about the Barbero-Immirzi parameter 7. In four- 
dimensional canonical loop quantum gravity, the Barbero-Immirzi parameter has been introduced 
in order to allow for the construction of the real (in the sense of the field of real numbers) su(2)- 
valued Ashtekar-Barbero connection \A7\ HH] . It is well known that it has no effects at the classical 
level (see |49H53] for a discussion about its controversial fate in the presence of fermions), since it 
drops out of the equations of motion by means of the Bianchi identities. However, in the quantum 
theory, it appears explicitly in the spectra of the geometric operators [SI [55] such as area and 
volume, and in the formula for the entropy of black holes |56H60j . It has been argued that the 
presence of this quantization ambiguity at the quantum level might be related to a "wrong" choice 
of the connection parametrizing the phase space at the classical level |61] . 
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In the context of four- dimensional spin foam models, the inclusion of the Barbero-Immirzi 
parameter has been made possible with the introduction of the EPRL and FK models. This has 
opened the possibility of relating the boundary states of spin foam models with the kinematical 
states of loop quantum gravity. However, several points still lack a clear explanation. In particular, 
the structure of the models seem to depend on the particular values taken by 7. For 7 < 1, the 
EPRL and FK models do coincide [22], and the space of boundary states can be reduced to the 
space of SU(2) spin networks. For 7 > 1 however, the two models are different, and even if the 
boundary states of the EPRL model coincide with the states of canonical loop quantum gravity, 
the space of boundary states of the FK model seems to be bigger. 

Since we have given in this work a proposal for a rigorous implementation of the spin foam 
quantization procedure in the three-dimensional context, we believe that it should be possible and 
interesting to investigate the issues arising in the presence of the Barbero-Immirzi parameter. In 
usual three-dimensional SU(2) Riemannian gravity, such an extension of the classical theory is 
not possible, since there is only a unique (invariant non-degenerate) bilinear form on su(2). For 
this reason, previous attempts to introduce a Barbero-Immirzi-like parameter were relying on the 
Chern-Simons formulation of BF theory with A 7^ 0, in which case there is another bilinear form 
available [65]. However, this does not mimic what happens in the four-dimensional theory, and 
renders difficult the comparison between the two situations. 



Since the three-dimensional Plebanski theory (1.4) is formulated for the gauge group S0(4), it 



can be generalized in a natural way to the case 7 7^ 0. Indeed, since there are two traces, denoted 
by Tr(- , •) and Tr(* • , •), available on so(4), it is possible to define two different pairings of elements 
of so (4), and to construct the modified action 

Sbf^[B,oj] = ^^xe^^'P (ty{B^,,F,p) + iTr(*S^, F.^)^ , 

where 7 G M — {0}. To see that 7 is the analogue of the four-dimensional Barbero-Immirzi 
parameter, and in particular that it disappears on-shell, let us split the B field into its self-dual 
and anti self-dual components, and write the action as 



where 



WSBF7Pi3,'*M= / d3xe'^'^^Tr((±)i3^,WF,p,. 



Following what we have done in section |lj it is possible to solve the torsion- free condition (if 

is invertible), and to plug back the solution for the connection in the action S'bf7 in order to obtain 

the second order Einstein-Hilbert action 



If we impose that B^^ be simple (i.e. if we choose the gravitational sector of solutions to the 



simplicity constraints), we know from (1.11a) and (1.12) that we have e"*" = e~ and ^^^g^u = ''~^9fj.u, 
and therefore it is clear the the Barbero-Immirzi parameter 7 disappears from the classical theory. 

This argument shows that it is possible to naturally extend the three-dimensional Plebanski 
action that we have introduced in order to have a non- vanishing Barbero-Immirzi parameter which 
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behaves classically exactly as in the four-dimensional theory. The Holst-Plebanski action 



Spi^[B, ooA] = \ I d^x L^'^'P f Tr(S^, F^p) + -Tr(*5^, F,p)] + c/)'"' Tr{B^,*B, 
^ Jm I \ 7 / 



can now be taken as a starting point for the spin foam quantization. It is clear that 7 plays no role at 
the classical level, and we can then study its fate in the quantum theory. It has been claimed in the 
four-dimensional context that a canonical quantization based on the non-commutative spacetime 
Lorentz-connection leads to a quantum theory independent of 7 |63) . It is therefore natural to ask 
if 7 can disappear as well in the spin foam quantization, and, if not, to investigate why it is not 
the case. We think that this question deserves a detailed treatment, and keep it for future work. 
We want to stress that this is the first time that a Barbero-Immirzi parameter is introduced in 
three-dimensional gravity. 



Discussion and conclusion 

In this work, we have introduced a three-dimensional Plebanski action for the gauge group SO (4), 
which is based on the notion of simplicity for one-forms. This simplicity condition on the B field is 
imposed via quadratic simplicity constraints, which admit, just like in the four-dimensional case, 
a gravitational, a topological, and a degenerate sector of solutions. We have shown at the classical 
level that the action reduces to the Einstein-Hilbert action in the gravitational sector, while in the 
topological sector the theory is trivial. The canonical analysis reveals that the theory is indeed 
equivalent to three-dimensional gravity. Interestingly, this model has a Hamiltonian structure 
similar to that of the four-dimensional Plebanski theory, and in particular features secondary 
second class constraints arising from the stabilization of the primary simplicity constraints. 

This three-dimensional Plebanski action allows to test the imposition of the simplicity con- 
straints in the spin foam quantization of gravity. We show that neither the BC nor the EPRL 
prescriptions lead to the expected vertex amplitude, which is that of the Ponzano-Regge model. 
However, the two schemes differ drastically in that only the EPRL model can give the proper 
amplitude when it is supplemented with the imposition of the secondary second class constraints. 
This suggests that the weak imposition of the simplicity constraints is an essential ingredient in 
the spin foam quantization, but that it is not sufficient in this three-dimensional model. We have 
given a concrete realization of the idea that imposing the secondary second class constraints is a 
missing element of the spin foam approach [33]. 

A limitation of the present model is that, because of its three-dimensional nature, it has trivial 
off-diagonal simplicity constraints, and therefore the BC and EPRL prescriptions lead to the same 
vertex amplitude (up to a normalization factor). This is a key difference with the four-dimensional 
case, where the off-diagonal simplicity constraints are non-trivial, and might very well encode the 
information about the secondary second class constraints. In order to investigate this possibility, it 
would be interesting to extend the present three-dimensional model to an arbitrary non-simplicial 
cellular decomposition. 

Interestingly, the three-dimensional SO (4) Plebanski action allows for the introduction of a 
Barbero-Immirzi parameter. This could be a good starting point to understand and test issues 
related to the role of the Barbero-Immirzi parameter in the canonical and spin foam quantizations 

m- 

Finally, let us conclude with a remark on the BC model. In the present work, we have followed 
only one very specific way to derive the BC model, which is the most straightforward one. There 
are however in the literature more rigorous and transparent approaches which lead to variants of 
the usual four-dimensional BC model, and which allow to overcome some of the criticism which 
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were originally raised against the model and to clarify its geometrical meaning. Details about these 
constructions can be found in the work of Bonzom and Livine [65] and Baratin and Oriti [31] . In this 
latter in particular, the authors use a non-commutative metric representation of group field theory 
to define a model of four-dimensional constrained BF theory, whose amplitudes reproduce a variant 
of that of the EC model. In this approach, it is shown that any spin foam model for constrained 
BF theory can be written as a simplicial path integral featuring an explicit non-trivial measure 
over the holonomies and the bivectors. The measure over the holonomies is naturally interpreted as 
imposing secondary constraints arising from the consistent imposition in each tetrahedron frame of 
the (primary) linear simplicity constraints. It would therefore be interesting to investigate wether 
these secondary constraints agree with those coming from the canonical analysis, and to test the 
proposal of j3T| [32] on the present toy model to see if it leads to the Ponzano-Regge amplitudes. 
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Appendix A: Notations 



In this work, notations are such that fj,,^, . . . refer to spacetime indices, a,b, . . . to spatial indices, 
I,J,... to so(4) indices, and to 5u(2) indices. We assume that the three-dimensional 

spacetime manifold A4 is topologically S x M, where S is a two-dimensional manifold without 
boundaries. We use the notation for vectors in with components (A*^, A*), the wedge product 
between elements u, u G (zero-forms or one-forms) to denote the operation {u A vY = e^j^^u^v'^, 
and the dot product for u • t> = u^Vi. We often denote the vectors A* S simply by A. 



Appendix B: The Lie algebra so (4) 

so(4) is the real Lie algebra of the isometry group S0(4) of the quadratic form 6 = diag(l, 1, 1, 1). 
We use capital Latin letters for internal vector indices and define the antisymmetric tensor e^-^^^ 
such that e*^^^^ = 1 and sijkl = s^"^^^- The indices are lowered and raised with the metric 6. 
The action of the Hodge dual operator is defined by 



, t/ J _ ^IJ jKL 



and it satisfies = id. 

One of the basis of so (4), composed of rotation generators Li and boost generators Xj, with 
i G {1,2,3}, has the following commutation relations: 

[Li, Lj] = Sij^Lk, [Ki,Kj] = Eij^Lk, [Ki, Lj] = Sij^Kk, 

where Sijk = ^^ijk- Starting from this basis of so (4), it is convenient to define a new basis **'Jj as 

^^^J^ = ^{L,±K,), 
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whose generators realize two commuting copies of su(2), i.e. satisfy 

For any element £ so (4), we have the decomposition = '+'^ + *~'^ into self-dual and anti self-dual 
components, where '^'^ = '*'^*'*'Jj. The action of the Hodge dual operator on the (anti) self-dual 
components if given by 

^{±)^ = ±{±)^. 

Finally, the vector representation of so (4) is given by 

Appendix C: Sign of det(*='='iJ) in the topological and gravitational sectors 

In this appendix, we compute the sign of det(**'i?) in the gravitational and topological sectors. 

1. Gravitational sector 

In the gravitational sector, we have Bj/ = s^'^xL^^ ^jl^ boost and rotational components 

of Bj/ are given respectively by 



2 

Therefore, we have the relation 

_ ,,-1 



and the self-dual and anti self-dual components ( l.lla| ) of the B field can be written as 



^^^Bl = Tix A e,y + ix'el - x'e',) = - = (-1 ^ Xo\)L',- (CI) 
Here we have introduced the three-dimensional matrix 

/ 




X 



associated to x such that xct* = (x A a)* for any a £ M'^, and I denotes the three-dimensional unit 
matrix. With this notation, we can compute from (|C1|) the determinant 



det((±)i?) = det(-I T Xo'x) det(4) = - (l + Xo^xl + xl + xD) det(L; 
Therefore, we see that in the gravitational sector we have e"*" = e^. 
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2. Topological sector 

In the topological sector, we have Bj/ = *£^'^klX^ ^jl = ^"^i ~ X"^ ^jj.^ and the boost and rotational 
components of Bj/ are given respectively by 

Kl = x%-X'el, V^ = {xf\e,)\ 
Therefore, we have the relation 

Ll = x,\x^K^)\ 



and the self-dual and anti self-dual components (1.11b) of the B field can be written as 

= -ix A e,y ± ix'el - x%) = " = " Xo'xW^. 
From this formula, we can compute the determinant 

det((±)i?) = det(Tl - xo'x) det{K;) = t(1 + Xo^xl + xl + xD) det(4). 
Therefore, we see that in the topological sector we have e+ = — e~. 

Appendix D: Algebra of constraints 

The algebra of constraints is 

{/Co, /Co} = 0, 

{/Co,/Ca} =Tr(^o,5a) ~0, 

{/Ca,/C{,} = 0, 

{/Co,^o} = -^o~0, 

{/Co,*a} =0, 

{/Ca,*o} = -^a + 2Tr(7ro,PaSo) « 0, 

{lCa,^b} =2TT{Tro,VaBh)^0, 

{/Co, /Co} = 0, 

{/Co,^a} =Tr(^o,*5a) ~0, 
{/Ca,-^o} = -TV(fi„,*7ro) «0, 

{/Ca,/Cfe} = 0, 

{/Co, <^ab} = -TT{VaBo,*Bb) - Tl{VbBQ,^Ba) = -^ab ~ 0, 

{/Cc, $ab} = -TT{VaBc,*Bb) - TT{VbBc,*Ba) + eacTr(7ro, [Bo,*Bb]) + e6cTr(^o, [Bo,*Ba]) « 0, 
{/C^,Coo} = -2Tr(S^,*5o) = -2C^o ~ 0, 

{iCf^Coa} = -Tl{B^,^Ba) = -C^a ~ 0, 
{/C/ijCaft} = 0, 
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^'o,*o} = 0, 

^'0,*a} = -2£^'^TV(Fcd,25a^0), 

^o,iCo} =£'''Ti{Fab,*Bo), 

^o,}Cc} = e^''Ti:{Fab,*Bc) - 2TV(Pc7ro, ^^o) « e^^ Ti:{Fab,*Bc), 
'^a,^C^,} ~0, 

*o,$afe} = -2Tr{VaBo,mBo) - 2TT{VbBo,*VaBo), 

*c,^afe} = -2Tr{VaBc,*VbBo) - 2Tt{VbB^,*VaBo) - e'^'Tv{Fde,eea[Bo,*Bb] + ecb[Bo,*Ba]), 
^'m'Coo} =0, 

^i„Coa} = 2Tr{i.B^,VaBo), 

^^^,Cab} = 2Tr{icB^,VaBb) + 2IV(*5^, PbS„). 
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